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Abstract. A novel symmetry method for finding exact solutions to nonlinear PDEs is 
illustrated by applying it to a semilinear reaction-diffusion equation in multi-dimensions. 
■ The method uses a separation ansatz to solve an equivalent first-order group foliation sys- 

C^hI tem whose independent and dependent variables respectively consist of the invariants and 

fH ' differential invariants of a given one-dimensional group of point symmetries for the reaction- 

diffusion equation. With this group-foliation reduction method, solutions of the reaction- 
diffusion equation are obtained in an explicit form, including group-invariant similarity 
solutions and travelling-wave solutions, as well as dynamically interesting solutions that are 
not invariant under any of the point symmetries admitted by this equation. 
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5 ! 1 Introduction 

The construction of group foliations using admitted point symmetry groups for nonlinear partial 
, differential equations (PDEs) is originally due to Lie and Vessiot and was revived in its modern 

r> ■ form by Ovsiannikov [T]. In general a group foliation converts a given nonlinear PDE into an 

. equivalent first-order PDE system, called the group-resolving equations, whose independent and 

dependent variables respectively consist of the invariants and differential invariants of a given 
one-dimensional group of point symmetry transformations. Each solution of the group-resolving 
equations geometrically corresponds to an explicit one-parameter family of exact solutions of 
the original nonlinear PDE, such that the family is closed under the given one-dimensional 
symmetry group acting in the solution space of the PDE. 

Because a group foliation contains all solutions of the given nonlinear PDE, ansatzes or 
differential-algebraic constraints must be used to reduce the group-resolving equations into an 
overdetermined system for the purpose of obtaining explicit solutions. Compared with classical 
symmetry reduction [21 [3], a main difficulty to-date has been how to find effective, systematic 
ansatzes that lead to useful reductions. 

An important step toward overcoming this difficulty has been taken in recent work [31[5] on 
finding exact solutions to semilinear wave equations and heat equations with power nonlineari- 
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ties. Specifically, this work demonstrates that the group-resolving equations for such nonlinear 
PDEs have solutions arising from a simple separation ansatz in terms of the group-invariant 
variables. Through this ansatz, many explicit solutions to the nonlinear PDE are easily found, 
whose form would not be readily obvious just by trying simple direct ansatzes using the original 
independent and dependent variables in the nonlinear PDE, or by simply writing down the form 
for classical group-invariant solutions. In particular, some of these solutions are not invariant 
under any of the point symmetries of the nonlinear PDE and thus fall completely outside of clas- 
sical symmetry reduction (while others coincide with explicit group-invariant solutions). Most 
importantly for applications, many of the solutions also have interesting analytical properties 
related to critical dynamics, blow-up behaviour, asymptotic behaviour and attractors. 

We will illustrate this group-foliation reduction method by applying it to obtain explicit exact 
solutions for the semilinear radial reaction-diffusion equation 

Ut = Urr + {n — l)r~^Ur + {p — ku'^)u, k = ±1, p = const (1) 

for u{t,r), with a nonlinearity power q ^ 0, —1, where r denotes the radial coordinate in n > 1 
dimensions or the half-line coordinate in n = 1 dimension. The symmetry structure of this 
reaction-diffusion equation is given by j6j 

time translation Xi = d/dt for all n, q, p, (2) 
scaling 'X.2 = 2td / dt + rd / dr — {2/ q)ud / du only for p = 0, (3) 

space translation X3 = d/dr only for n = 1, (4) 

where X is the infinitesimal generator of a one-parameter group of point transformations acting 
on {t,r,u). For constructing a group foliation, it is natural to use the time translation genera- 
tor ([2]), since this is the only point symmetry admitted for all cases of the parameters n > 1, 
q 7^ 0,p. 

In Section [21 we first set up the system of group-resolving equations given by the time- 
translation symmetry ([2]) for the reaction-diffusion equation ([T]), which uses the invariants and 
differential invariants of the symmetry generator Xi as the independent and dependent variables 
in the system. We next state the form required for solutions of the group-resolving system to 
correspond to group-invariant solutions of the reaction-diffusion equation ([1]) with respect to 
the point symmetries generated by Xi, X2, X3. 

In Section [3l we explain the separation ansatz for directly reducing the system of group- 
resolving equations. This reduction yields an overdetermined system of differential- algebraic 
equations which can be readily solved by computer algebra. We present the explicit solutions 
of these equations and then we derive the resulting exact solutions of the reaction-diffusion 
equation. These solutions include explicit similarity solutions in the case p = 0, and explicit 
travelling wave solutions in addition to an explicit non-invariant solution in the case re = 1. 

In Section IH we show how the success of the reduction ansatz can be understood equivalently 
as constructing partially-invariant subspaces for a nonlinear operator that arises in a natural 
way from the structure of the group-resolving equations. This important observation puts our 
method on a wider mathematical foundation within the general theory of invariant subspaces 
developed by Galaktionov [7]. 

Finally, we make some general concluding remarks in Section [5j 

2 Group-resolving equations and symmetries 

To proceed with setting up the time-translation group foliation for the reaction-diffusion equa- 
tion ([1]), we first write down the invariants (in terms of t, r, u) 



X = r, V = u, 



(5) 
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satisfying Xix = Xif = 0, and the differential invariants (in terms of ut, Ur) 

G = Ut, H = Ur, (6) 

satisfying ^^^G = X^^^i/ = where X^^'' is the first-order prolongation of the generator ([2]). 
Here x and v are mutually independent, while G and H are related by equality of mixed r, t 
derivatives on ut and Ur, which gives 

DrG = DtH, (7) 

where Dr, Dt denote total derivatives with respect to r, t. In addition, v, G, H are related 
through the reaction-diffusion equation ([1]) by 

G - r^~'^Dr{r'^~^H) = {p - kv'i)v. (8) 

Now we put G = G{x, v), H = H{x, v) into equations ([7]) and ([8]) and use equation ([5]) combined 
with the chain rule to arrive at a first-order PDE system 

G^ + HG„ - GH^ = 0, (9) 
G-{n- l)H/x - H^- HHy = {p - kvi)v, (10) 

with independent variables x, v, and dependent variables G, H. These PDEs are called the 
time-translation- group resolving system for the reaction-diffusion equation ([1]). 

The respective solution spaces of equation ([T|) and system ([9]), (fTO|) are related by a group- 
invariant mapping that is defined through the invariants ([5]) and differential invariants ([6]). 

Lemma 1. Solutions {G{x,v), H{x,v)) of the time-translation- group resolving system Q, (|10|) 
are in one-to-one correspondence with one-parameter families of solutions u{t, r, c) of the reac- 
tion-diffusion equation ([T]) satisfying the translation-invariance property 

u{t + e, r, c) = u{t, r, c(e, c)), 
where c(0, c) = c in terms of an arbitrary constant c and parameter e, such that 

ut = G{r,u), Ur = H{r,u) (11) 
constitutes a consistent pair of parametric first-order ODEs whose integration constant is c. 

We now examine the relationship between the symmetry structure of the reaction-diffusion 
equation ([TJ and the symmetry structure inherited by the time-translation-group resolving sys- 
tem dHD, (doD. 

Firstly, through the identifications defined by the variables ([5]), (l6|), the prolongation of any 
point symmetry generator X = oiXi + 02X2 + 03X3 of equation ([1]) has a natural projection 
to a point symmetry generator Y = 02X3^^ + 03X3^^ modulo x[^^ of system ([9]), (llOp . The 
time-translation Xi thus gets annihilated by this projection, i.e. Yi = 0, while the scaling X2 
and the space-translation X3 respectively project to 

Y2 = xd/dx - {2/q)vd/dv - 2(1 + l/q)Gd/dG - (1 + 2/q)Hd/dH when p = 0, (12) 

and 

Y3 = d/dx when n = 1. (13) 
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Secondly, with respect to these inherited symmetries (fT2|) and p^ . the system ([9]), (fTO|) has 
a reduction to ODEs yielding solutions where (G, H) is invariant respectively under scalings 

X — )• Arr, V — )• \~'^l'^v when p = 0, 
and under translations 

X — >• X + e when n = 1. 
Thus, translation-invariant solutions have the form 

{G,H) = {g{v\h{v)) (14) 
satisfying the ODE system 

(/i/5)' = 0, g-hh' = {p-kv'i)v. (15) 
Scaling-invariant solutions have the form 

(G, i?) = (x-2-2/'?(^(y), x-^-^Z-^/iC^^)), V = vx^l'i (16) 
satisfying the ODE system 

((/i + 2y/g)/5)' = -2^, (7 + (2 - n + 2/g)/i - (/i + 2y/g)/i' = -A;y^+^ (17) 

Integration of the parametric ODEs (fTT]l for such solutions (fT4|l , (fT5]l and (fT6]l , (fT7|l leads to the 
following two correspondence results. 

Lemma 2. In the case n = 1, there is a one-to-one correspondence between solutions of the 
translation- group resolving system Q, (|10p with the invariant form (|14|) and one-parameter 
families of travelling-wave solutions of the reaction- diffusion equation ([T|) given by the group- 
invariant form u = /(^) where, modulo time-translations t ^ t -\- c, the variable ^ = r — t/a 
is an invariant of the translation symmetry X = ad/dt + d/dr = aXi + X3 in terms of some 
constant a [determined by ODE system (jlSp ). 

Lemma 3. In the case p = 0, there is a one-to-one correspondence between solutions of the 
translation- group resolving system ([9]), (jlOp with the invariant form (jl6p and one-parameter 
families of similarity solutions of the reaction- diffusion equation ([T]) given by the group-invariant 
form u = r~'^/'^f{^) where, modulo time-translations t ^ t + c, the variable ^ = t/r"^ is an 
invariant of the scaling symmetry X = 2td/dt + rd/dr — [2/q)ud/du = X2. 

Furthermore, in all cases, static solutions u{r) of the reaction-diffusion equation correspond to 
solutions of the translation-group resolving system with G = 0. Hereafter we will be interested 
only in solutions with G 7^ 0, corresponding to dynamical solutions of the reaction-diffusion 
equation. 

3 Main results 

To find explicit solutions of the group foliation system Q, (jlOp for {G{x,v), H{x,v)), we will 
make use of the same general homogeneity features utilized in [31[5]. First, the non-derivative 
terms {p — ku'^)u in the reaction-diffusion equation ([TJ appear only as an inhomogeneous term 
in equation ([TO]) . Second, in both equations ([9|) and (fTOj) the linear terms involve no derivatives 
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with respect to v. Third, the nonlinear terms in the homogeneous equation Q have the skew- 
symmetric form HG^ — GH^, while HH^ is the only nonlinear term appearing in the non- 
homogeneous equation ([TO]) . Based on these features, this system can be expected to have 
solutions given by the separable power form 

G = gi{x)v + g2ix)v°' , H = hi{x)v + h2{x)v"' , a / 1. (18) 

For such a separation ansatz (|18p . the linear terms G^, G, H/x^ in equations @ and (|10|) 
will contain the same powers u , v"" that appear in both G and and moreover the nonlinear term 
HGi, — HyG in the homogeneous equation ([9]) will produce only the power v"" due to the identities 
v'^{v)y — {v'^)vV = (a— and v{v)v — {v)yV = v'^{v'^)y — {v^)yv'^ = 0. Thus, equation ([9]) can be 
satisfied by having the coefficients of v and v'^ separately vanish. Similarly the nonlinear term 
HHy in the non-homogeneous equation (jlOp will only yield the powers v, t)", u^""^. Since we 
have a 7^ 1 and q ^ 0, equation (llOp can be satisfied by again having the coefficients of v and v'^ 
separately vanish and by also having the term containing 7;2a-i ij^lance the inhomogeneous 
term kv''^^ . In this fashion we find that equations ([9]) and (llOp reduce to an overdetermined 
system of 5 algebraic-differential equations for gi{x)^ 52(2^)) hi{x), h2{x), together with the 
relation a = 1 + q/2. This system can be solved by a systematic integrability analysis. We have 
carried out this analysis using the computer algebra program Crack [8], which contains a wide 
repertoire of modules for reduction of differential order and polynomial degree, splittings with 
respect to the independent variable(s), eliminations, substitutions, factorizations, integrations, 
and length-shortening of equations, among others. 

Proposition 1. For q ^ 0, —1, the separation ansatz (llSp yields altogether five solutions of the 
translation- group resolving system ([9]), (jlOp with G 7^ 0; 

G = ±(4 - n) (j^'^) ^("-')/("-'V^, 

H = ± (j^fc) yin~mn-2) + (2 - n)v/x 

p = 0, n/2,3, g = 2/(2-n); (19) 
G = ^{pv + {pk)'/'v'^'^/'), 

H = ±((-^Y\+(-^]' v'^'^/'], n = l, g/-2; (20) 
G = 

H = ±l{^) v-{^) v'^'^/'], n = l, q^-2; (21) 

G = ^{pv± {pkf^ tanh {{p/2)^l^{x + ci))v'^), 

H = ±{k/2y/^v^ + (p/2)i/2 tanh {{p/2y/\x + ci))v, 

n = l, q = 2, p>0; (22) 
G=^{pvT{-pk)^/Han{{-p/2)^/\x + ci))v^), 
H={± ik/2)^/\^ - (-p/2)i/2 tan {i-p/2)^/\x + ci))v), 

n = l, q = 2, p<0. (23) 

Solution (|19p satisfies the scaling-invariance reduction l\l<j\) : solutions l\20\) and l\21\) satisfy the 
translation-invariance reduction (IE 
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Remark 1. A shift of the arbitrary constant ci — t- ci it i7r/(2p)^/^ in solution (|22p and ci 
ci lb 7r/(— 2p)^/^ in solution (|23|) respectively yields the two further solutions 

G=^{pv± {pk)^/^ coth ((p/2)V2(x + ci))v^) , 

H = ±{k/2y/^v^ + (p/2)i/2 coth {{p/2y/^{x + ci))v, 
n = 1, 9 = 2, P > 0, 

and 

G=^{pv± (-pA:) cot ((-p/2)i/2(x + ci))^;^), 

H={± ik/2)^/\^ + (-p/2)i/2 cot {i-p/2)^/\x + ci))t;), 
n = 1, 9 = 2, P < 0, 

as obtained through the trigonometric identities 

tanh(6l ± i7r/2) = coth(6'), tan(6' ± 7r/2) = - cot (61). 



The form of the separation ansatz (jlSp can be naturally generalized to include additional 
powers 

N N 

G = gi{x)v + Y,9i+i{xK', H = hiix)v + Y,hw{x)v''\ (24) 

i=l 1=1 

where 

Oj / I, Oi - Oj / 0. (25) 

This multi-term ansatz (for > 2) leads to a considerably more complicated analysis compared 
to the previous two-term ansatz (with = 1). Specifically, the homogeneous equation Q now 
contains the powers y'^i+'^j-'^ in addition to v, v'^\ while the non-homogeneous equation (jlOp 
contains the further powers v'^'^^~^, and v'^'^^. To determine the exponents in these powers, 
a systematic examination of the possible balances (which rapidly increase in number with A^) is 
necessary. We carry out this analysis by computer algebra, again using the program Crack [8]. 
In particular, for any fixed A > 1, Crack can be run automatically by setting up a priority 
list of modules to be tried recursively with case splittings given a high priority, in order to split 
equations Q and (jlOp under the ansatz (j24p . ()25p and in each case solve the resulting overde- 
termined system of algebraic-differential equations for gi and hi, i = 1,2, . . . , N . Importantly, 
the modules in Crack are able to organize the case distinctions in an efficient way that lessens 
a combinatorial explosion by using algebraic conditions (equalities and inequalities) arising in 
the steps of solving the equations to keep track of whether a possible balance of two powers of v 
generates a new case (or subcase) which has to be solved or not. Running Crack, we find that 
the system dlD, ([10]), ([25]) has solutions when A^ = 2 but not when N = 3. (Based on this 
outcome we have not tried to investigate the system when A > 3.) 

Proposition 2. For q ^ 0,-1, and A < 3, the separation ansatz (j24p yields two additional 
solutions of the translation- group resolving system (jlOp with G ^ 0: 

G = ?,{pv + {k'pf/-' + {kp'f'\^l^), 

H = -{2/pY/\pv + {epY/^ + {kp')'/\'/'), 

n = l, 9 = -3/2, p>0; 
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G = A{pv + {-pk)^/^v^^/^ + (-VA:)V4t;i/3), 

H = -{3/p)'/^pv + {-pk)'/^v-'/' + i-ip'k)'/\'/'), 
n = l, q = -8/3, p > 0. 

Both solutions satisfy the translation-invariance reduction (jl4p . 

We now obtain explicit solutions u{t, r) of the reaction-diffusion equation ([T]) from the solu- 
tions {G{x, v), H{x, v)) of its translation-group resolving system ([9]), (fTOj) by integrating the cor- 
responding pair of parametric first-order ODEs (llip . This integration yields a one-parameter so- 
lution family u{t, r, c) which is closed under the action of the group of time-translations t — )■ t + e. 

Theorem 1. The semilinear reaction- diffusion equation ([T]) has the following exact solutions 
arising from the explicit solutions of its translation- group resolving system found in Proposi- 
tions 1 and 2 (and Remark 1); 

y V(^-3)(n-2)y V2 r 
n^2,3, g = 2/(2-n), p = 0; (26) 



u ^ i- ik/pf'^ + exp Uq (^) ± (, + 4) (^) it + c) 



n = l, q^-2, p^O; (27) 
u = {{k/pfl' + exp ( Tg ^7^77^1 k ± (9 + 4) f t^tT^I + c) ) ) | , 



2(g + 2) J \ ^ ' V2(g + 2), 
n = l, (/^-2, (28) 
n = /(r-3(p/2)^/2(t + c)) where 

In (p/(e) + {kp^f'^'fiif^ + {k'pf'^) - (2/^/3) arctan((2(p/A;)i/V(0'/' + 1)/^^) 

= -(2p)i/2^, n = l, g = -3/2, p>0; (29) 
n = /(r-4(p/3)^/2(t + c)) where 

In [pfiif^ + (-pfc)i/2 + {-Ap'^kf/^fiif'^) - 2arctan ((-4p/A;)i/V(6^^^ + l) 

= -4(p/3)i/2e, n = l, g = -8/3, p>0; (30) 
cosh((p/2)V2(r + ci)) 



u 



{k/py/^ sinh((p/2)V2(r + ci)) ± exp(-3p(t + c)/2) ' 
n = l, g = 2, p>0; (31) 
sinh((p/2)V2(r + ci)) 



u 



{k/p^/^ cosh((p/2)i/2(r + a)) ± exp(-3p(t + c)/2) ' 
n = l, g = 2, p>0; (32) 
cos((-p/2)V2(r + ci)) 



u 



-/c/p)i/2 sin((-p/2)V2(r + ci)) ± exp(-3p(i + c)/2) ' 
n = 1, g = 2, p < 0; (33) 
sm{{-p/2y/'^{r + ci)) 



u 



-k/pyl'^ cos((-p/2)V2(r + ci)) =F exp(-3p(t + c)/2) ' 
n = 1, g = 2, p < 0, (34) 

where c, ci are arbitrary constants. 
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Remark 2. Solutions (j3ip and (j32p are related by shifts of the arbitrary constants ci — t- ci zt 
m/{2pY/'^ and c — t- c=F ^7r/(3p) through the trigonometric identities 

sinh(6'±i7r/2) = ±«cosh(6'), cosh(6l±OT/2) = ±zsinh(6'), 

while solutions ()33p and ()34p are related in a similar way through the trigonometric identities 

sin(6'±7r/2) = ±cos(6l), cos(6l±7r/2) = =Fsin(6l). 

Modulo time-translations, solution (]26p is a similarity solution, since it has the form shown in 
Lemma El This solution has been obtained in our previous work [5] using a two-term separation 
ansatz to solve the group-resolving equations given by the scaling symmetry ^ for the reaction- 
diffusion equation ([T]) in the case p = 0. 

Solutions (j27p to ()30p are travelling- wave solutions with the form shown in Lemma [2] (where 
the constant a is given by the wave speed). The two solutions ()27p and ()28p were first obtained 
in through standard symmetry reduction. To the best knowledge of the authors, the other 
solutions ([29]) and pOj) are new. 

In contrast, solutions (j3ip to ()34p are non-invariant solutions, as they contain both t and r 
whereas the symmetry group of the reaction-diffusion equation ([T]) in the case p 7^ and n = 1 
is generated by both time-translation and space-translation symmetries ^ and These 
solutions have been found previously in [U [10] using Bluman and Cole's nonclassical method. 



4 Invariance properties of the reduction ansatz 



As we will now explain, the success of the separation ansatzes (llSp and (124p for reducing the 
system of group-foliation equations Q and (jlOp has a mathematically natural interpretation 
within the theory of invariant subspaces [7]. 

We start by writing the group-foliation equations in an evolutionary operator form 

h) ^ = '^{{11)) = {{kv'i - p))) + (g + (1 - n)H/)) + ( -HH, ") ' ^^^^ 

where ^ defines a nonlinear operator acting on the pair of variables (G, H). Then we view the 
ansatz (llSp as defining a linear space of functions with a separable power form 



where the coefficients depend only on x. This linear space (j36p is not invariant under <I> because 

\^^{ (a-l)(Ml-M2) \ a 

51 + (1 - n)hi/x -h\- p) \g2 + (1 - n)h2/x - {a+ l)/ii/i2/ 

produces terms that have additional powers v*^^^. 

However, the operator $ will preserve a set of functions contained in the linear space (|36p if 
the power a and the coefficients gi, g2, hi, /i2 satisfy the conditions 



2a — 1 = q + 1, a/i2 = k, 
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so that the terms containing v and v'^'^ cancel out. In this sense the linear space will then 
be partially-invariant [7] under <I>. As a consequence of this invariance property, the evolution 
equation ([35]) for (G, H) reduces to an overdetermined system of differential-algebraic equations 
which we can solve for gi, g2, hi, as discussed in Section [2j 

The separable power ansatz (|18|) thus can be viewed as explicitly constructing a partially- 
invariant linear subspace for the group- foliation operator <l>. A similar discussion applies to the 
multi-term ansatz ()24p . 

It is worth pointing out that apart from the mathematical interpretation given to the an- 
satzes (jlSp and (j24p . the theory of (partially) invariant linear subspaces does not provide any 
general constructive method or algorithm for finding a successful ansatz. So the main aspect of 
the present work is to demonstrate the effectiveness of a general separation ansatz for solving 
group-foliation equations. We also emphasize that this approach is much more effective than if 
an analogous separation method were to be applied either directly for solving the given reaction- 
diffusion PDE or instead for seeking a (partially) invariant linear subspace for that PDE itself. 

5 Concluding remarks 

In general the method of group foliation reduction using a separation ansatz as illustrated in this 
paper is able to yield exact solutions to nonlinear 2nd order PDEs with power nonlinearities. 

This method works with any admitted group of point (or contact) symmetries and gives 
a systematic reduction of the group foliation equations into an overdetermined system consisting 
of algebraic equations and 1st order differential equations that can be derived and in most cases 
solved by means of computer algebra (e.g. using the program Crack [8]). In particular, for 
a given nonlinear 2nd order PDE having two independent variables, solutions are produced in 
an explicit form, whereas standard symmetry reduction only gives a 2nd order ODE that still 
has to be solved to find group-invariant solutions explicitly and in general this step can be quite 
difficult. 

Moreover, because the group foliation equations contain all solutions of the given nonlinear 
PDE, the method can yield solutions that are not invariant under any of the point (or contact) 
symmetries admitted by the nonlinear PDE. 

It is straightforward to extend this method to higher-order PDEs with power nonlinearities. 
As well, it should be possible to apply the same method to PDEs having more general forms 
of nonlinearities by utilizing a general separation of variables ansatz with respect to all of the 
independent variables in the group foliation equations. In general, the success of such an ansatz 
can be interpreted mathematically [7] as constructing a (partially) invariant linear subspace for 
a nonlinear operator coming from the structure of the group foliation equations. 

Related work using a similar group-foliation method applied to nonlinear diffusion equations 
appears in [12] . Group foliation equations were first used successfully in [131 El US] for obtaining 
exact solutions to nonlinear PDEs by a different method that is applicable when the group of 
point symmetries of a given PDE is infinite-dimensional, compared to the example of a finite- 
dimensional symmetry group considered both in [H [5] and in the present work. 
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